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We theoretically propose an experimental method to probe electronic excitations in graphene - 
monoatomic layer of carbon - by monitoring the fluorescence quenching of a semiconductor quantum 
dot (or a dye molecule) due to the resonance energy transfer to the graphene sheet. We show how 
the dispersion relation of plasmons in graphene can be accurately extracted by varying the back-gate 
voltage and the distance between the quantum dot and graphene. 



I. INTRODUCTION 

Many appealing properties of graphene - a 
monoatomic crystalline sheet of carbon - stem from its 
unique electronic structure.^ Specifically, its honeycomb 
lattice combined with the conjugation of 7r-electrons 
over the entire sheet results in the electronic spectrum 
of a zero-gap semiconductor with "ultra-relativistic" 
electrons and holes .^^I Already this makes graphene 
enormously appealing from both basic and applica- 
tion standpoints. What makes graphene even more 
attractive is the possibility to tune these properties in 
a wide range by patterning, chemical functionalization, 
doping etc. For example, shifting the Fermi level away 
from the charge neutrality (Dirac) point by applying 
the back-gate voltage, and, therefore, changing the 
graphene's electrical cond uctiv ity, can become useful in 
graphene-based electronics .'i'^ 

The improved electrical conductivity of a back-gated 
or chemically doped graphene sample leads to qualitative 
changes in its optical properties. In particular, collective 
excitations (plasmons), rather then single-particle exci- 
tations (electron- hole pairs), are expected to define an 
electronic response of graphene to a low-frequency op- 
tical perturbation. Plasmons in graphene are of great 
interest from the basic perspective, as collective excita- 
tions in a two-dimensional electron gas with very pecu- 
liar properties. Besides, the graphene plasmonics holds 
promise for, e.g., pho toni c and optoelectronic devices, as 
well as metamaterials.l^Plasmons in graphene have been 
predicted and studied, theoreticallySHlIl^j^^ experimental 
studies of this phenomenon are still very sparse.'^^'^ 

In this paper we propose an experimental technique 
to probe plasmons in graphene. The technique is based 
on the Forster resonance energy transfer between a fluo- 
rescent semiconductor quantum dot (or a dye molecule) 
and a nearby graphene layer. Efficient energy trans- 
fer between a dye molecule and graphene has been the- 
oretically predicted by Swathi and SebastiarP^Wl^ and 
later confirmed experimentally.^^ However, Swathi and 
Sebastian described electronic excitations in graphene 



on the single-particle level, which is accurate only for 
a nearly undoped (charge-neutral) graphene. In this pa- 
per, electronic excitations are treated more accurately 
by adopting the random-phase approximation, which al- 
lows for recovering the collective electronic behavior. We 
demonstrate that electronic excitations in graphene, both 
single-particle and collective, can be sensitively probed 
by studying the fluorescence quenching of the quantum 
dot due to the energy transfer to graphene. Specifically, 
we show how the plasmon dispersion can be extracted 
from experiment. 

The expected advantage of the proposed technique 
over th e typi cally used electron energy loss spectroscopy 
(EELS)i3HISlig its intrinsic locality, i.e., plasmons are pro- 
posed to be probed locally by a semiconductor quantum 
dot (typically a few nanometers in diameter). In con- 
trast, EELS averages the plasmonic response over a cer- 
tain portion of a graphene sample, thus, adding the in- 
homogeneous broadening due to, e.g., charge puddles,!^ 
to experimental observables. 

The paper is organized as follows. The general theory 
of quantum dot fluorescence quenching due to the reso- 
nance energy transfer to graphene is given in Sec.|ll] The 
analysis of quenching efficiency within the single-particle 
and random-phase approximation levels is provided in 
Sec. IV and Sec. Ill respectively. Sec. [V] concludes. 



II. 



FORSTER RESONANCE ENERGY 
TRANSFER 



Forster resonance energy transfer (FRET) refers to 
the transfer of electronic excitation energy between 
chromophore^^U rnediated by the nonradiative Coulomb 
coupling.l^^l This process consists of deexcitation of an 
initially (optically) excited donor chromophore and the 
simultaneous excitation of an acceptor chromophore. An 
example of such a process is the energy transfer between a 
semiconductor quantum dot (QD) and a non-fluorescent 
organic molecule as an acceptor chromophore.^'^ As the 
energy transfer to the non-fluorescent "dark" molecule 
(i.e., quencher) competes with the intrinsic fluorescence 
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of the "bright" QD, the FRET rate can be assessed 
through the effective decrease (quenching) of the QD flu- 
orescence quantum yield, and/or through the shorter ap- 
parent fluorescence lifetime, measured by time-resolved 
fluorescence spectroscopy. Specifically, the apparent flu- 
orescence lifetime is given by r = 1/ {tq^ + fcg), where tq 
is the lifetime of the isolated QD and kq is the FRET 
(quenching) rate. Accordingly, the decreased fluores- 
cence quantum yield in the presence of FRET is given 
by t/tq. Therefore, kg can be extracted provided t and 
To are known from experiment. 

The efficiency of FRET is strongly affected by param- 
eters of the acceptor's excitation spectrum. A trivial ex- 
ample is a vanishing FRET rate in a system where an 
acceptor chromophore does not have excitations resonant 
to the lowest excited state of a donor chromophore. This 
allows one to use FRET as a spectroscopic tool to probe 
an electronic structure of a system of interest. An advan- 
tage of this method is the near-field regime of the donor- 
acceptor interaction, which allows for probing excitations 
forbidden in the far-field (optical) regime due to, e.g., the 
large wavelength mismatch between optical photons and 
material electronic excitations. Short-wavelength charac- 
ter of electronic excitations in graphene, especially that 
of plasmons,!^^ hinders the study of these excitations by 
standard far-field optical techniques. We propose to use 
FRET between QD and a graphene sheet to probe elec- 
tronic excitations in graphene. The theoretical frame- 
work for FRET in the QD-graphene complex is developed 
in the rest of the section. 



A. Fluorescence quenching efficiency 

We define the wavefunctions of the excited and the 
ground states of QD as |e) and \g), respectively. The true 
many-body ground state and excited states of graphene 
are denoted by |0) and |n), respectively. The quenching 
rate, kq, for the energy transfer from QD to graphene is 
given by Fermi's golden rule as 



2nh-'J2\(n\{g\V\e}m'S{e-E„ 
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where V is the operator of Coulomb interaction between 
fluctuating charge densities of QD and graphene 



V= / drV{r){\e){g\ + \g){e\)p{r)- 



(2) 



The operator of graphene charge density is given by 
— e/5(r) — —e(f\r)ip{r), where operators (^^(r) and (p{r) 
create and destroy an electron at position r within the 
graphene sheet, respectively. The absolute value of the 
electron charge is denoted by e = |e|. The excitation 
energies of QD and graphene are denoted by e and En, 
respectively. Vector variables are denoted in bold. 

For the rest of the paper, we adopt the dipole approx- 
imation for QD, which results in V{r) = — e(d • r)/r^. 



The transition dipole of QD is given by d. This approx- 
imation is accurate if z 3> D, where z is the distance 
between QD and graphene, and D is the QD diameter. 
The validity of this approximation for the realistic case 
of PbSe QD is discussed in Sec. [IV] 

The quenching rate in Eq. ([I]) can be rewritten through 
the retarded polarization operator of graphene IV {q, e) as 
(a detailed derivation is given in the Appendix) 
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where dy and d±^ are the projections of the QD transition 
dipole d onto the graphene's plane and the normal to this 
plane, respectively. If we ignore for simplicity that the 
population of the QD excited state depends on the angle 
between the QD transition dipole and the polarization 
of the initial QD-exciting laser pulse, we simply need 
to average over all the possible orientations of d with 
respect to the graphene's plane which yields (c?|)/2 = 

{d\) = resulting in 



kq — 



qdq Im [11'' (q, e)] 
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(4) 



We define the quenching efficiency (QE) as ipq — kg/ kr , 
where kr is the fluorescence rate for the isolated QJJ^ 



4eJ 



(5) 



Substituting Eq. ([5| into Eq. Q, one obtains for QE 



qdq Im [TI''{q, e)] e 
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(6) 



B. Massless Dirac fermions approximation 

Polarization operator n'"(g, e) for graphene can be ob- 
tained at various levels of theory. At zeroth order approx- 
imation with respect to the electron-electron Coulomb 
interaction within the graphene sheet, it can be evalu- 
ated adopting the free massless Dirac fermions (MDF) 
approximation. At this level, the polarization opera- 
tor, denoted by flj^q, e), is a bare polarization bubble 
describing a single non-interacting electron-hole pair in 
graphene. Within this approximation, IIJKq, e) has been 
evaluated previouslj^^"^ for arbitrary doping level and its 
imaginary part is shown in Fig. [l](a) . This figure is valid 
for any level of doping, defined by the chemical poten- 
tial, /i, measured relative to the Dirac point, since due 
to the linear dispersion relation of massless fermions in 
graphene, the polarization operator does not change with 
doping, if measured in units of and plotted against 

unitless coordinates q/kf and e/n, where Vf « 0.4 is the 
Fermi velocity in atomic units and kf = ^/hvf^^ 

With respect to the wave number, q, the MDF approx- 
imation is only valid when g <C a^^, where a sa 2.46A is 
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FIG. 1: Density plot of (a) Im[n5(g,e)] and (b) 
Im [n^P^(g, e)] in units of y^^. Checkerboard pattern em- 

phasizes regions where the imaginary part of the polarization 
operator vanishes exactly. 



the lattice constant of graphene. Because of the factor 
g-2gz ^jjg integrand of Eq. this condition is equiv- 
alent to restricting z ^ a, i.e., to considering the quench- 
ing efficiency ipq only at large QD-graphene distances 
compared to the lattice constant of graphene. However, 
because of the adopted dipole approximation we are al- 
ready restricted to z ^ Z?, where the typical QD diame- 
ter {D) is on the order of a few nanometers. Thus, setting 
z ^ a because of the MDF approximation does not re- 
strict the range of applicabihty of Eq. ^ any further 
compared to the one set by the dipole approximation. 
With respect to the excitation energy, the MDF approx- 
imation has been shown to be valid for e < 1 eV.^^ 

The key features of Im [nQ(g, e)] include (i) the singu- 
larity along the e/^ = q/kf line, which corresponds to 
single-particle excitations with (q, e) vector lying within 
the surface of the Dirac cone, and (ii) the absence of 
single-particle excitations, i.e., Im[no(g, e)] = 0, in re- 
gions A and B, marked by checkerboard patterns in 
Fig. l|^a). Equation 4] with ^^^{q, e) substituted in is very 
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fluorescence quenching of a dye molecule due to single- 
particle excitations in graphene. 

The polarization operator within the bare bubble ap- 
proximation does not include the graphene's polarization 
self-consistently which can become crucial at non-zero 
doping levels (/i > 0), where the finite carrier density at 
the Fermi level leads to the efficient Coulomb screening 
within the graphene's sheet. To correct for this, we eval- 
uate the polarization operator within the random-phase 
approximation (RPA) a^^ 



n 



RPA 



l~W{q)ni{q,e) 



(7) 



2i:e /kq is the two-dimensional 



where W{q) 

Fourier transform of the Coulomb potential within the 
graphene's plane. The effective dielectric constant of the 
environment is denoted by h. A free-standing graphene 
sheet in vacuum corresponds to k = 1. For graphene 
laying on top of a half-space dielectric substrate with 
dielectric constant k, the effective constant is given by 
K = (k + l)/2.^^^^ For example, for a Si02 substrate 
(k=4) the effective dielectric constant is k = 2.5. In the 
rest of the paper, except for results shown in Fig. [4j the 
vacuum conditions (k = 1) are assumed. 

The imaginary part of n^p^((7, e) is depicted in 
Fig. [T|b). It is seen that the "single-particle" singularity 
e//i — q/kf is gone and the new singularity at e oc q^^"^ 
appears instead. This emergent singularity corresponds 
to the collective electronic excitation, i.e., plasmon, in 
graphene.'SHU It is rather pedagogical to see exactly how 
this singularity appears from Eq. ^j. The naive sub- 
stitution of nQ((7, e), shown in Fig.^l] into Eq. ([t]) re- 
sults in Im [njjp^(5, e)] vanishing exactly in the entire 
region A since both nQ(g, e) and W{q) are real in this re- 
gion. However, the more careful analysis reveals that the 
imaginary part of Hq(5, e) is not exactly zero in A, but 
infinitesimal instead due to the usual small imaginary 
constant in the denominator of the Lindhard function, 
which preserves causality. The infinitesimal imaginary 
part of the bare polarization operator can be safely ne- 
glected if HQ(g, e) is analyzed by itself. However, when 
1 — M^(q)Re [HQ(g, e)] vanishes in the denominator of 
Eq. I?]), this small imaginary part yields singularity (i5 
function) in Im [Hpp^(g, e)], i.e.. 



Im [n5jp^(g, e)] 5{q - qp{e)), 



(8) 



t of Swathi and Sebastian'i^ for the 



where the plasmon dispersion relation is given by (j'p(e) ~ 

at small e, i.e., at e//i ^ 1, and obtained numerically 

at higher excitation energies .'^Therefore, the general pre- 
scription to assume all the infinitesimal parameters finite 
till the very end of calculations proves to be critical in 
this case. 

It is important to note that since there is no single- 
particle excitations in region A, the plasmon lifetime is 
infinite [5 function in the imaginary part of Hpp^(g,e)] 
within this region.!^ Once plasmon "leaves" region A 
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FIG. 2: Regimes of asymptotic dependence of quenching ef- 
ficiency on the QD-graphene distance, at z — > oo (Regimes 
I, II and III) and at 2: — 5> (Regime IV) . 



{q/kf > 0.65), it acquires the finite lifetime due to the 
Landau damping and soon disappears. 

In the following two sections, |TTT]and|lVj we discuss the 
efficiency of the QD fluorescence quenching due to un- 
screened excitations in graphene, described by nQ(g, e), 
as well as screened excitations in graphene, described by 
n^jPA (<?'€), respectively. 



III. QUENCHING BY UNSCREENED 
EXCITATIONS 



QE due to unscreened single-particle excitations in 
graphene can be obtained from Eq. Q by substituting 
W{q,e) with no(g, e). This approximation is accurate if 
either (i) the effective dielectric constant of a substrate 
is high, and, therefore, the effective electron-electron 
Coulomb interaction within graphene is greatly reduced, 
or (ii) e ^ /Lt, i.e., graphene becomes effectively undoped. 
In both cases, the screening within the graphene's layer 
becomes ineffective yielding n5jp^(g, e) ~ ^o{q, e)- 

Four different regimes of the asymptotic dependence of 
QE on the Ql^graphene distance, z, are shown schemat- 
ically in Fig. [2] At large distances, exponent e"^"^^ in 
the integrand of Eq. ^ decays rapidly with q, which 
guarantees that at fixed excitation energy e the domi- 
nant contribution to QE comes from lowest possible g's 
where the imaginary part of the polarization operator is 
still nonzero. The two qualitatively different cases are (i) 
"gapless" , where the imaginary part of the polarization 
operator is already non-zero even at infinitesimally small 
q > 0, and (ii) "finite-gap", where Im[ng(g, e)] becomes 
non-zero only at some finite q, i.e., at q > q* with finite 



q* > 0. The typical power law dependence of ^^{q, e) on 
q at g — ?► in the gapless case 



and 



at q 



and 



nji(g,e)cxg", 

n5(g,e)(x(g~g*)", 
g* + in the finite-gap case, yield 
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at z — >■ cx), respectively. In what follows, we consider 
three different asymptotic regimes of QE at z — >■ 00 
(Regimes I, II and III). The asymptotic behavior of QE 
at z — > (Regime IV) is described subsequently. 

Regime I. At high excitation energies (e//i > 2) the 
imaginary part of IlfKg, e) is seen in Fig. IJa) to corre- 
spond to the gapless case. Further, it can be shown that 
Im[nQ(q, e)] cx q^ at small q (see Eq. (12) in Ref. 7 for 
the long wavelength limit of the polarization operator), 
which is then combined with Eq. (10a) to give 



iPq CX 1/z^. 



(11) 



This power-law dependence is shown as Regime I in 
Fig. [2] The asymptotic behavior in Eq. (11) is satis- 



fied at e//z > 2, and, therefore, also at e /i, where 
^ can be treated as zero, i.e., the graphene sheet be- 
comes effectively undoped. The asymptotics of l/z"* for 
the fluorescence quenching by single-particle excitations 
in undoped (and also weakly dope d) gr aphene was first 
predicted by Swathi and Sebastian! ^^*^''! 

Generally, the asymptotic dependence 1/z^ is rather 
expected since QE due to the energy transfer from a chro- 
mopliore to a quencher accompanied by the interband ex- 
citation of single electron-hole pairs in an iV-dimensional 
quencher typically scales as z^^^, where z is the distance 
between the chromophore and the quencher. FRET to 
Oc? quencher (e.g., small organic molecule) corresponds 
to TV = 0, naturally resulting in 1/ z^W^ Nanowire as a 
quencher {N = 1) leads to 1/z^ asymptotics.!^ Excita- 
tion of single-electron holes in a dipole-to-surface con- 
figuration, i.e., when a quencher occupies the half-space 
(A^ = 3), yields 1/z'^ asymptotics.!^ The graphene sheet 
is a 2d object, which naturally leads to 1/z^ dependence 
of QE at large z. 

Regime II. At 1 < e//i < 2, the imaginary part of 
nQ(q, e) vanishes exactly at g < g*, where the threshold 
value oi q* /kf = 2 — e//z marks the onset of interband 
single-particle excitations. It can be shown that in this 
finite-gap case, Im [nQ(g, e)] is proportional to {q^- g 
at g — >■ g* -|- 0, which yields [by virtue of Eq. ( 10b )] 
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which is marked as regime II in Fig. [2j 



5 



Regime III. At low excitation energies, e//i < 1, the 
finite-gap case is again realized with q* /kf = e//i, which 
is the onset of intraband single-particle excitations. In 
this regime, Im [nQ((j', e)] oc (g — resulting in 



-l/2-2q'z 



(13) 



depicted as regime III in Fig.[2| 

Regime IV. Finally, at a fixed excitation energy there 
are no single-particle excitations with q/kf > 2 + e/fi, i.e, 
in region B in Fig. [T] This introduces the natural high-g 
cutoff for integration in Eq. If, at certain (small) 

z, exponent e^^'^^ is still w 1 at g approaching this cut- 
off, then QE becomes constant with respect to z. This 
is depicted as regime IV in Fig. [2j However, we expect 
this regime to be hardly accessible experimentally. For 
example, at e = 0.8 eV and /i = 0, the realization of this 
regime requires z <5 A. At such small distances photoin- 
duced charge transfer from QD to graphene can become 
the dominant mechanism of fluorescence quenching pro- 
hibiting the analysis of the relatively less efficient energy 
transfer channel. Besides, at such distances both dipole 
and MDF approximations are likely to break down. 
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FIG. 3: (a) Dependence of the quenching efficiency (p^ on 
the QD-graphene distance z. (b) The filled area shows where 
ipq{z) is between 1000 and 1 at a given excitation energy e//i. 
The red (left) and black (right) circles mark ipq{z)=1000 and 
ipq{z)=l contour lines, respectively. 



IV. QUENCHING BY SCREENED 
EXCITATIONS 



To evaluate QE in the case of screened excitations in 
graphene, one has to substitute W{q,e) in Eq. ^ with 
n^p^((7, e) given by Eq. ([t]). First, we consider regime 
I in Fig. [2j In this regime, the product W{q)IlQ{q,e) 
becomes proportional to 1/qx q^ — q at q ~> 0, which re- 
sults in the approximate equality n^p^(g, e) w IlQ{q,e) 
at small q. Thus, the asymptotic behavior of >fq{z) at 
large z is the same for screened and unscreened excita- 
tions in regime I; i.e., taking screening into account does 
not lead to qualitative changes in QE. This can be eas- 
ily understood for e ^ /i, where graphene becomes ef- 
fectively undoped, and, therefore, the small free carrier 
density renders screening within graphene inefficient. 

The situation is different in regime III (e//i < 1), where 
taking Coulomb screening into account within RPA leads 
to the emergence of the new (plasmon) singularity in the 
polarization operator, as is shown in Fig. p|b), region 
A. Equation (jsj) implies that Im [n5jp^(g, ej] vanishes 
exactly at q < q* — qp{e)\ i.e., the finite-gap situation is 
realized with q* defined by the plasmon dispersion rela- 
tion. Substituting Eq. ([8| into Eq. ([6]), one obtains 



(14) 



The absence of the power-law multiplier in front of the 
exponent, which was universally present in the finite-gap 
situations in the previous section, is related to a delta- 
functional instead of a power-law singularity of the po- 
larization operator. 

The asymptotic behavior, given by Eq. ( |14[ ), is cor- 
rect even outside regime III, since the plasmon branch in 



Fig. [T](b) remains singular up to e//i ~ 1.3, i.e., well 
within what used to be regime II in the case of un- 
screened excitations. For 1.3 < e//i < 2 the imaginary 
part of Wj^p^{q,e) scales as {q - g*)^/^ at q ^ q* + 
with the q gap defined hy q* /kf — 2 — e/ii, which yields 
ipq{z) OC z~^/^e~^' ^, i.e., the large-z asymptotics is iden- 
tical to that of regime II in the case of unscreened exci- 
tations. 

To examine how accurately large- 2: asymptotics for 
^Pq{z) reproduce exact solutions at finite z, we numeri- 
cally evaluate the integral in Eq. ^ for the realistic case 
of PbSe QD with the excitation energy of e = 0.8 eV. 
At this energy, the fluorescence quantum yield of the iso- 
lated PbSe QD can be as high as 70-90% with fluores- 
cence lifetimes up to 1 ^s,'^"' suggesting that QE can be 
directly extracted from experiment, since the observable 
fluorescence quantum yield is given by 1/(1 -f (^g) in the 
presence of FRET.^e 

Numerically evaluated ^q{z) for several values of chem- 
ical potential in the range /i = 0.2 — 1.6 eV is shown 
in Fig. 3]^a). The smallest value of chemical potential 
adopted corresponds to Regime I since e/ fi — A. The cor- 
responding ^Pq(z), depicted by black circles, is expected 
to show l/z** dependence at large z, and, indeed, demon- 
strates slowly decaying non-exponential tail. 

All the other values of chemical potential give e//i < 2, 
and, therefore, are expected to give ^Pq{z) oc e"^'''''^'^^^ 
and ^Pq{z) cx z~^/^e~^'' ^ at large z for e/// < 1.3 and 
1.3 < e//i < 2, respectively, with q* defined by q* /kf = 
2 — e/^. As expected, all (/?q(z)'s, except for the one 
corresponding to the lowest value of chemical potential 
(/i = 0.2 eV), demonstrate the nearly exponential decay 
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FIG. 4: QE decay rate q* plotted vs e/^. Circles (black), 
squares (red) and triangles (blue) correspond to single- 
particle approximation, RPA for free-standing graphene (vac- 
uum) and RPA for graphene on substrate (Si02, k=4), re- 
spectively. 



at large z. The filled area in Fig. [sj^b) shows the range 
of QD-graphene distances, where QE is between 1 and 
1000 - somewhat loosely chosen range where the accurate 
experimental measurement of QE is still possible. QE is 
seen to decay the fastest with z (lowest z at fixed ipq) 
where the q-gap is the largest, i.e., at e//i w 1.3, as is 
seen in Fig. ^h). 

The rate of the exponential decay of QE at e//i < 2 
directly reflects the specific structure of Im [n^p^(g, e)] 
at small q, i.e., the width of the finite g-gap. There- 
fore, the large-z behavior of QE can be used to extract 
the valuable information about electronic excitations in 
graphene, e.g., the plasmon dispersion. To illustrate 
how the plasmon dispersion can be extracted, we plot 
q* against e/^ in Fig. |4j where q* is the decay rate of 
QE at large z, obtained by fitting the large-^ QE decay 
with (fq{z) oc e"^'' ^. Specifically, q* is plotted for the 
case of unscreened excitations (black circles), as well as 
for screened excitations in free-standing graphene (red 
squares) and graphene laying on top of the Si02 sub- 
strate (blue triangles). Scanning through the range of 
values of chemical potential using the back-gate, we are 
able to extract the dispersion relation of singularities of 
the imaginary part of the polarization operator. As is 
seen, the linear dispersion relation of single-particle ex- 
citations (black circles at e//j, < 1) as well as that of the 
plasmon in the free-standing graphene (red squares at 
e/^ < 1.3) in graphene is accurately recovered. 

As discussed in the beginning of Sec. electronic 
excitations in graphene are effectively single-particle ei- 



ther if (i) there is a strong substrate-induced dielectric 
screening or (ii) the excitation energy is high, i.e., e ^ /z. 
The first case is illustrated by the plasmon dispersion for 
graphene on the Si02 substrate (blue triangles), which 
is closer to the single-particle excitations (black circles) 
than that corresponding to the free-standing graphene 
(red squares). The second case is effectively realized at 
e/n > 1.3, where dispersions of electronic excitations 
with and without accounting for the in-graphene screen- 
ing are nearly identical. 



The proposed method to extract the plasmon dis- 
persion requires an accurate control over the distance 
between PbSe QD and the graphene sheet. This can 
be accomplished by either using the core-shell type-1 
structures'^! or by growing a dielectric layer on top of 
graphene by atomic layer deposition with the controllable 
thickness and then depositing QDs on top of that dielec- 
tric layer. The first approach can be based on PbSe/CdSe 
core-shell structures, where the large bulk bandgap of 
CdSe (-1.7 eV) as compared to that of bulk PbSe (~0.28 
eV) guarantees the confinement of both electron and hole 
within the PbSe core (type-1 structure) This confine- 
ment guarantees that the shell serves only as an inert 
spacer between the PbSe core and the graphene layer. 
Recent advances in core-shell structure fabrication tech- 
niques allow one to control the thickness of the shell with 
the monolayer (subnanometer) precision.l^S Atomic layer 
deposition provides an alternative strategy for control- 
ling the distance betw;een QD and graphene with sub- 
nanometer resolution. 



Finally, we wish to discuss the effect of the finite plas- 
mon propagation length on the applicability of the pro- 
posed method. So far, we assumed that the plasmon has 
infinite propagation length within region A in Fig.[l|^b). 
However, additional damping channels, not accounted 
for in RPA (e.g., defect scattering, electron-phonon cou- 
pling), can "broaden" the plasmon in region A. It fol- 
lows from Eq. (6|, that if the plasmon width in the q 
domain, Sq, is less then 1/z, then Im [n^p^(g, e)] can 
still be treated as 6 function, resulting in Eq. (14). Ac- 



cordingly, deviations from the exponential dependence 
are expected at z > Sq^^ ~ I, where I is the plasmon 
propagation length. 



Our numerical tests show that the exponential decay 
of QE is typically already established at z comparable 
with the plasmon wavelength A. On the other hand, es- 
timations by various authors suggest that the plasmon 
propagation length could be as high as 10-100 A! ^^ I '^° I 
Therefore, we expect QE to decay exponentially at z > A, 
thus allowing for extraction of plasmon dispersion. The 
onset of deviations from this exponential decay at larger 
QD-graphene distances, z > I, naturally provides a good 
estimate for the plasmon propagation length. 
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V. CONCLUSION 

Based on the detailed analysis of fluorescence quench- 
ing efficiency in the QD-graphene complex, we have pro- 
posed a method of probing and studying electronic exci- 
tations in graphene. The method has been demonstrated 
to be sufficiently sensitive to allow the extraction of the 
dispersion relation of plasmon in graphene. We hope that 
this study will stimulate experimental efforts in this di- 
rection, especially because the proposed method is based 
on the QD-graphene complex which can be of interest 
not only as a means to probe electronic excitations in 
graphene, but also on its own merit as a key component 
of hybrid nanostructures with promising properties. For 
example, the ability to excite plasmon locally in graphene 
using a semiconductor quantum dot can become of great 
use in graphene plasmonics. Another recently proposed 
use of QD-graphene complexes is in photovoltaics.^ 

The RPA-based description of electronic excitations 
in graphene, adopted in this paper, does not include 
such many-body effects as exchange and correlation, im- 
purity and defect scattering, or electron-phonon cou- 
pling. These effects can affect the plasmonic response 
of graphene in two ways. First, there can be deviations 
of the actual plasmon dispersion relation from the one 
shown in Fig. [T]^b). Second, additional channels of plas- 
mon damping can appear, as discussed in the previous 
section. The proposed experimental technique is capable 
of assessing both the dispersion relation and the damp- 
ing and, thus, is expected to stimulate the development 
of beyond-RPA theoretical methods by providing a nec- 
essary experimental validation. 
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Appendix: Derivation of quenching rate 

Evaluating explicitly the QD part of the matrix ele- 
ment in Eq. ([T]) one obtains 



kg^2TTh-^J2 ("I / drV{r)p{r)\0} 

J a 



S{e-E,,). (A.l) 



This can be rewritten as 

^9 = (I^E / ^^^^'/ dqdq'F*(q)e-'i"-y(q')e*'>''-' 

(A.2) 



(27r)3 



where Pno(r) = (n|/5(r)|0) and V(r) = 
/ c?q y(q)e*'i ''. Further, the 5 function can be sub- 
stituted using the identity (5(e-£;„) ^ ^ J die*('-^")*/'^ 
yielding 



kq = I dvdv' I dqdq'V*{q)e-^^''V{q')e''^ 



/..n>(r,r';.)e-^ 



(A.3) 



-iE„t/h 



IS 



where n>(r, r';t) = -iH ^ I]„ Pon(r)/Ono(r')e 
the "greater" polarization operator for graphene in the 
Lehmann representation. At sufficiently large distances 
between the quantum dot and the graphene layer, ^(r) 
varies smoothly within the graphene's plane. This makes 
it possible to average the polarization operator over the 
unit cell with respect to both r and r'. After this av- 
eraging, the polarization operator becomes insensitive 
to variations of electronic density on the scale of the 
graphene's unit cell and, therefore, acquires the isotropy 
and the continuous translational symmetry instead of 
the discrete one, leading to a possibility to substitute 
n>(r, r';t) n>(|r — r'\,t). Then, integrations over r, 
r', and t can be interpreted as spatial and time Fourier 
transforms, respectively, resulting in 



ih 
J2 



^ |dq|y(q)l'n>(g,e), (A. 



4) 



where 



n>{q,e)^ J dr J dtn>(|r|,i)e~*('i'"-<^*/''). (A.5) 

Using the relations between real-time correlation and re- 
sponse functions at equilibrium,^^ which at zero temper- 
ature yields 11^ = 2ilm[n^], we obtain 



f^i^-J^ J dci\V{q)flm[U^iq,e)]. (A.6) 



This equation is similar to the well-known result by Metiu 
[Eq. (2.30) in Ref.gS]. 

Finally, we provide without derivation the two- 
dimensional Fourier transform (within the graphene's 
plane) of V{r) 

F(q) = 27rie(d|| cos(6i) -|- ^d_L)e"'^ (A.7) 

where d|| is the projection of the QD transition dipole 
onto the graphene's plane, d\\ — |d|| | being its magnitude. 
The angle between vectors dy and q is denoted by d. The 
projection of the QD transition dipole onto the normal 
to the graphene's plane is denoted by d±. Substitution of 
V{q) into Eq. (A.6 1 and the subsequent integration over 
e yields Eq. 
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